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Abstract 

Aggregation processes with an arbitrary number of conserved quantities are investigated. 
On the mean-field level, an exact solution for the size distribution is obtained. The asymp- 
totic form of this solution exhibits nontrivial "double" scaling. While processes with one 
conserved quantity are governed by a single scale, processes with multiple conservation 
laws exhibit an additional diffusion-like scale. The theory is applied to ballistic aggre- 
gation with mass and momentum conserving collisions and to diffusive aggregation with 
multiple species. 
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I. Introduction 



Irreversible aggregation processes under ly many natural phenomena including, e. g., 
polymerization [1], gelation [2], island growth [3-4] and aerosols [5]. The classical rate 
theory of Smoluchowski describes the kinetics of such processes [6-9]. Recently, scaling 
[10-16] and exact [17-23] theoretical studies showed that spatial correlations play a crucial 
role in low dimensions. While the above examples are diffusive driven, there are physical 
situations such as formation of large-scale structure of the universe [24] and clustering 
in traffic flows [25], where the aggregates move ballistically. So far, theories of ballistic 
aggregation [26] have been restricted to scaling arguments [26-30]. 

In the ballistic aggregation process both the mass and the momentum are conserved. 
In polymerization processes involving copolymers, each monomer species mass is a con- 
served quantity. Hence, we study aggregation processes with multiple conservation laws. 
The simplest example for such a system is aggregation with k distinct species. Both the 
multivariate distribution and the single variable distributions are of interest. We present 
exact solutions to the time dependent and the steady state mean-field rate equations. Al- 
though they are straightforward generalizations to the well known solutions they exhibit 
interesting behaviors. An asymptotic analysis shows that fluctuations associated with a 
single conserved quantity are Gaussian in nature. As a result, an additional "diffusive" 
size scale emerges. 

We apply the above theory to ballistic aggregation as well as diffusive aggregation. In 
the case of ballistic aggregation, we use an approximate collision rate to obtain a solution to 
the Boltzmann equation in arbitrary dimension. While this approach agrees with the scal- 
ing arguments, it suggests that for a given mass, the momentum distribution is Gaussian. 
We compare these predictions with one and two dimensional simulations. Furthermore, 
we consider steady state properties of the aggregation process by introducing input of 
particles. For homogeneous input, a novel time scale describing the density relaxation is 
found. In the case of a localized input, clustering occurs only for d < 2. Additionally, 
we apply the theory to two-species aggregation with diffusing particles. Using the density 
dependent reaction rate, we obtain the leading scaling behavior of the two relevant mass 
scales. 

This paper is organized as follows. In section II, we present exact solutions of the 
rate equation theory. We investigate time dependent as well as steady state properties of 
the process. We then apply the theory to ballistic aggregation with momentum conserving 
collisions (section III) and to diffusive driven aggregation (section IV). We conclude with 
a discussion and suggestions for further research in section V. 
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II. Theory 

Following the above discussion, there are aggregation processes where several physical 
quantities are conserved. In such processes, it is natural to label the aggregates by a "mass" 
vector m = (mi, • • • , rrik), where every component represents a conservation law. Let us 
denote the probability distribution function for particles of mass m at time t by P(m,t). 
Mean field theory of the binary reaction process assumes that reaction proceeds with a rate 
proportional to the product of the reactants densities. Thus the mean field approximation 
neglects spatial correlations and therefore typically holds in dimensions larger than some 
critical dimension d c [10, 14-16]. The rate equations [6] are written as 

^P(m',t)P(m-m',t)-2P(m,t)^P(m',t), (1) 

where the sum is carried over all k variables m^, 1 < i < k. The loss term prefactor reflects 
the fact that two particles are lost in each collision. One can verify that the rate equations 
conserve each mass separately, %. e., ^mjP(m, t) = ^miPo( m ) 5 where Po(m) is the 
initial distribution. In writing Eq. (1) we have implicitly assumed that the rate if(m, n) 
at which the reaction (m) + (n) — > (m + n) proceeds does not depend on masses of the 
reactants, if (m, n) = const. This constant can be set to unity without loss of generality. 
For an arbitrary reaction rate kernel K(m,n), Eq. (1) is easily generalized so that, e. g., 
the gain term is replaced by Yl K(m! ', m — m')P(m / , t)P(m — m', t). 

To solve Eq. (1) we introduce the generating function, F(z,£), defined by 

F(z,t)= 5> m P(m,t). (2) 

{mi} 

In the above equation we have used the shorthand notations z = (zi, ■ • • , Zk) and z m = 
z™ 1 ■ ■ ■ z™ h . Moments of the distribution function are readily obtained by evaluating the 
generating function and its derivatives in the vicinity of the point z = 1 = (1, . . . , 1). For 
example, the total cluster density is given by N(t) = F(l,£). The equation describing 
the temporal evolution of the generating function F(z,t) can be evaluated by a proper 
summation of the rate equation. This yields 

— — = F 2 — 2FN. (3) 

dt ' 

As a preliminary step, we evaluate the time dependence of the density. The corresponding 
rate equation is obtained by evaluating Eq. (3) at z = 1, 

IE = ~ N - (4) 
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Without loss of generality we set the initial density to unity and therefore we have N(t) = 
1/(1 + t). We also note that by subtracting Eq. (4) from Eq. (3), a simple differential 
equation for the quantity F — N, emerges, d(F — N)/dt = (F — N) 2 . This equation is 
readily solved to find 

^■'>= (1+t ) ( ffrU W) - (s) 

where the notation Fq(z) = F(z,t = 0) has been used. Eq. (5) represents the general 
solution to Eq. (1). Indeed, it is a simple generalization of the well known solution [6] for 
the single mass aggregation process. 

We consider the simplest multivariate case, 

k 

Po(m) =k- 1 ^6(m i -l)l[6(m j ). (6) 

i=l jjti 

These initial conditions imply F = z = (z\ + h Zk)/k and from Eq. (5) we have 

F(z,t) = - r-T- -. (7) 

V ; (l + t)(l + t-tz) K 1 

A possible application of these initial conditions is to aggregation processes involving k 
distinct species with equal initial densities. Although this symmetric situation appears 
to be too simple at first, it contains the necessary ingredients for exploring the long time 
kinetics. 

Before we investigate the multivariate distribution, let us first study properties of the 
following distribution function 

P(m+, t)= P ( m ' ^ 6 ( m + -(rn 1 + --- + m k )). (8) 

{mi} 

This distribution corresponds to the sum variable m + = mi + ■ ■ ■ + rrik- It is useful to 
introduce the generating function F(z, t) = z m+ P(m +1 1). This generating function can 
be obtained from F{z,t) by replacing z with z. Expansion of F(z,t) = z/(l + t)(l+t — tz) 
in powers of z yields the sum distribution 

fm + -l 

The variable m_|_ ignores the "identity" of the different conserved variables and thus, the 
problem reduces to ordinary aggregation. In the long time limit, m_|_ ~ t, and the corre- 
sponding distribution is given by 

P(m+,t) ~ t" 2 exp(-m+/t). (10) 
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Using the scaling variable, M + = m + /t, this distribution can be also written in terms of 
a scaling function P(m + ,t) ~ t~ 2 ^{M + ), with ^{x) = exp(— x). 

The multivariate probability distribution function can be now found by expanding 
F(z,t) and comparing with Eq. (2), 

P(m, t) = P(m+, t)g(m) with g(m) = g^>±l[ . (n) 

Note that ^| m .^(m) = 1. The above expression is the explicit solution to the mean field 
equations. However, its long time nature is of particular interest and we proceed with an 
asymptotic analysis. 

To study the asymptotic properties of P(m, t), we concentrate on the case k = 2 and 
then generalize the results to arbitrary k. The time independent geometric factor reads, 
^(mi, = 2 _m + (mi + m2)!/mi!m2! ~ m + 1//2 exp ( — (mi — m2) 2 /2m + ) . The limit of 
large masses, mi, 777.2 ^> 1, is the relevant one since the masses grow indefinitely. Using 
the limiting form of P(m + , t), we arrive at the asymptotic form of the mass density, 

P(mi,m2,t) ~ £ _2 m + exp(— m+/t) exp(— m 2 _/2m + ), (12) 

with m± = mi ± m-2. Two mass scales govern the mass distribution, m+ ~ t and |m_| ~ 
y/m+ ~ \ft. Furthermore, introduction of the scaling variables, M + = m+/t and M_ = 
m- /y/i, enables us to write the solution in a convenient scaling form 

P(m 1 ,m 2 ,t) ~ t" 5 / 2 $(M+,M_). (13) 

The scaling function $ depends on the scaling variables only, 

$(x, y) = x~ 1/2 exp(-x) exp(-y 2 /2x). (14) 

Interestingly, each of the variables, mi and m.2, exhibit simple scaling: mi,rri2 ~ t. 
Additionally, the moments M n (t) of the distribution function, M n (t) = X]{ mi } mn -P( m 5 1), 
are described asymptotically by linear exponents. By taking derivatives of the generating 
function, one can show that M n (t) ~ t a ^ with ct(n) = n + — 1 . This behavior agrees with 
ordinary scaling and therefore one could naively expect single-size scaling to hold. However, 
from the complete form of the mass distribution we learn that it is impossible to write the 
scaling solution in terms of a simple scaling function, P(mi, 777,2) ~ t~ 2 Q(mi/t,m2/t). 
Such a scaling solution would imply that the process has only one intrinsic size scale. 
Instead, this problem exhibits "double-scaling" , as the process is governed by two distinct 
size scales. The second scale |m_| ~ \ft is hidden, e. g., it does not clearly appear in 
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the moments of the distribution function. The mass difference is reminiscent of a diffusive 
scale for the following reason. The mass difference, m_ = mi — m 2 , is also a conserved 
variable, since each of its components is conserved. Hence, for an aggregate of total mass 
m + , the mass difference is a sum of m + random variables. According to the central limit 
theorem this variable is Gaussian and thus, m_ ~ y/m+- 

The above results can be generalized to arbitrary k. The generalized mass difference 
variable is m 2 _ = k~ x £\ (m^ — rrij) 2 and the mass distribution follows the scaling form 

P(m, t) ~ t- (fc+3)/2 $(M+, M_), (15) 

with the scaling function 

3>(x, y) = x~ (fc_1)/2 exp(-x) exp(-y 2 /2x). (16) 

The above results can be also generalized to situations with asymmetric initial con- 
ditions, i. e., the initial conditions are not invariant with respect to a permutation of the 
mass variables, {mj. Denoting by Mi = X]{m i } m *-^o(m) the i th conserved mass, one can 
easily show that with the transformation m$ — > rrii/Mi and P(m) — > Mi • • • MfcP(m), the 
system reduces to the symmetric case. For this transformation to be valid, the process 
must be truly multivariate, Mi 7^ 0. In this case, the prefactor g(m) equals probabilities 
associated with a biased random walk in a /c-dimensional "mass" space. 

Previous results have been established for a constant rate kernel, if (m, n) = const. 
Although the most general situation is hardly tractable, in the case of sum-variable depen- 
dent reaction rate, if (m, n) = K(m+, n+), the sum- variable distribution function P(m+, t) 
satisfies a simpler single-variable Smoluchowski equation. In this case, the complete multi- 
variate distribution function P(m, t) is related to P(m+, t) by Eq. (11). Therefore solvable 
variants of the single-variable Smoluchowski equation provide exact solutions for mul- 
tivariate Smoluchowski equation with the corresponding reaction rates. For the usual 
Smoluchowski equation, exact solutions have been found for three types of reaction rates: 
K(m, n) = const, K(m, n) = m + n, K(m, n) = mn, and for their linear combinations 
[6,5,31]. In particular, for the sum-kernel K (m, n) = |(m+ +n + ), subject to the monodis- 
perse initial conditions of Eq. (6), we again arrive at the exact solution of Eq. (11) with 

rm+fl-e-^l™*" 1 
P(m + ,t)=exp[-t-m+(l-e- t )] L +l (w+) , JJ . (17) 

Similarly, one can find an exact solution for the product kernel. 

We now investigate the steady state properties of the aggregation process by intro- 
ducing input of particles into the system. For simplicity, we consider homogeneous input 



6 



with rate h and restrict ourselves to the constant rate kernel. The governing equations are 
modified by adding an input term, 

dP ^ Q = p ( m '> t) p ( m - m '> *) - 2p ( m > *) p ( m '' *) + ^( m )- ( l8 ) 

The input function, P(m), satisfies the normalization condition Yl{m }-^( m ) = 1 smce the 
total input /i^{ mi }P( m ) is equal to h. 

The solution to the above equation parallels the solution to the time dependent prob- 
lem. We denote the steady state distribution by P(m) and the steady state generating 
function by P(z). This generating function is obtained by eliminating the time derivative 
of Eq. (18) , F 2 - 2FN + hR = 0, where R = R(z) = ^ {mi} z m P(m). The normalization 
condition R(l) = 1 is satisfied by the input function, and consequently the solution to the 
generating function reads 

F(z) = Vh (l - y/l - Rjzfj . (19) 

The total monomer density is given by = F(l) and thus, 

N(h) = Vh. (20) 

In analogy with the initial conditions of Eq. (6), we consider the input R(m) = Po(m), 
which in turn implies R(z) = z. In this case, the steady state generating function is simply 
F(z) = \fh(l — \/l — z) . Expanding in powers of the variables Zi, and comparing with the 
definition of the generating function, one finds the steady state distribution 

P(m) = P(m+Wm), with P(m + ) = ( 2m +V- ^1) 

v 1 v +;yv ; ' v +; (2m+ - l)(2 m +m+!) 2 V 1 

There is a strong similarity with the time dependent counterpart, P(m, t). The distribution 
among aggregates with the same total mass, m+, is given by the same combinatorial factor 
g(m) which has already appeared in solution to the time dependent problem. For large 
masses, the distribution describing the sum variable, m + , is an algebraic one, N~ 1 P(m + ) ~ 

3/2 

m + . We can write now the asymptotic form of the complete steady state distribution, 

P(m) = P(m+, m_) ~ Vh m+ {k+2)/2 exp(-m 2 _/2m+). (22) 

For the time dependent problem, an additional size scale emerges for the mass difference 
variable, m_ ~ y/rn+- In the steady state, on the other hand, the total mass diverges and 
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thus, the central limit theorem does not apply. However, the similarity between the two 
cases is still strong as for a fixed m+, the mass difference distribution function is Gaussian. 

For completeness, we briefly discuss the relaxation properties of the density towards 
the steady state. It is possible to obtain these properties by incorporating the time depen- 
dent density, N ~ t~ a , and the steady state solution, N ~ h? , into a single expression. 
Let us assume the scaling form [17] 

N(h,t) ~t" a ^(t/r), (23) 

with the relaxation time r ~ . In the limit of a vanishing input rate the time dependent 
solution must be recovered and thus, ip(x) — > 1 as x — > 0. In the long time limit steady 
state is approached and thus, if)(x) ~ x a for x 1. Comparing with the definition of the 
relaxation time the exponent relation 

7 = a(3 (24) 

is found. The exponent (3 = 1/2, is obtained from the decay exponent, a = 1, and the 
steady state exponent, 7 = 1/2. Hence, the relaxation time diverges in the limit of a van- 
ishing input rate, r ~ /i -1 / 2 . This result can be obtained directly by solving the equation 
dN/dt = —N 2 + h which is readily performed to find N(h, t) = v^tanh ((* + t )Vh). The 
time shift to is determined by the initial density, N(t = 0) = \^htanh(y^hto). We prefer 
the above scaling argument since it is applicable to a wide class of problems. 
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III. Ballistic Aggregation 

Ballistic aggregation is a natural and important example of a process with multiple 
conservation laws. Both the mass and the momentum are conserved quantities and thus, 
there are d + 1 conservation laws in d dimension [26] . A collision between two particles 
results in an aggregate whose mass as well as momentum are given by a sum over its 
components. One can view this system as a gas of sticky particles, i. e., an inelastic gas with 
a vanishing restitution coefficient. Heuristic arguments predict certain scaling properties 
of the system. However, other aspects of this problem, such as the mass-momentum 
distribution function remain unsolved. Our theory is well suited for this problem and 
yields an approximate form for the particle distribution function. 

We start by discussing the one dimensional case and then generalize to higher di- 
mensions. We denote the probability distribution function for particles of mass m and 
momentum p at time t by P(m,p) (we suppress the time variable). The stosszahlansatz 
Boltzmann equation, which describes the temporal evolution of this density, must conserve 
mass and momentum. Additionally, the collision rate between two particles is given by 
their velocity difference and hence, the Boltzmann equation reads 

- — ~^r~^ = \v' — v"\ P(m' ,p')P(m — m' ,p — p') — 2P(m,p) \v — v'\ P(m',p'), 

m' ,p' m' ,p' 

(25) 

where v' — v" = p' jm! — (p — p') /(m — m') and v — v' = p/m —p'/m'. This approximation 
ignores possible spatial correlations and is generally uncontrolled. By replacing the kernel 
terms \v— v'\ and \v'—v"\ with their average value (v), an approximate Boltzmann equation 
is written, 

dP ^ p ^ = (v) (J2 p ( m '' p') p ( m ~ m 'iP- p') - 2p ( m > p) p{m '' p')) • ( 26 ) 

The time-dependent factor (v) can be absorbed into a novel time variable, T, defined by 
dT/dt = (v). The resulting Boltzmann equation is equivalent to Eq. (1). 

For simplicity we consider an initial system of identical particles with zero average 
momentum. Without loss of generality we assume that typical initial quantities such as 
the mass and the momentum equal unity. Hence, the initial distribution function, P (m,p), 
is given by 

P (m,p) =6(m-l)(6(p-l) + 6(p+l))/2. (27) 

In principle, the general solution of Eq. (26) can be obtained using the generating function 
method. However, with these specific initial conditions, the process reduces to the two-mass 
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process discussed in the preceding section. Indeed, by identifying m with the total mass 
m + = mi +TO2 and p with the mass difference m_ = mi — 777-2, the initial conditions Eq. (6) 
and Eq. (27) are the same. Moreover, identical rate equations describe the time evolution 
of both P(mi, 777,2) and P(m,p). From the solution of Eq. (11), the mass-momentum 
distribution function is found in terms of the time variable T, 

Tm_1 2" m m! 
to m ) ~ ( 1 + T )m+i (( m+p )/2)!((m-p)/2)!' ( - 1 

Using the asymptotic scaling properties of the solution 777 ~ T and p ~ y/T, we can 
rewrite the solution in terms of the observable time t. Since dT/dt ~ u ~ p/m ~ T -1 / 2 , 
one has T ~ t 2 / 3 . Hence, the well known the scaling laws 

m~t 2/3 and p ~ t 1/3 , (29) 

are recovered. Asymptotically, the mass-momentum distribution is given by the following 
form 

P(p,m) ~ (m)- 2 m- 1 / 2 exp(- 7 ^- - |-Y (30) 

V (m) 2m/ 

with the average mass (m) ~ t 2 / 3 . As discussed in the previous section, the momentum 
p is a sum of m independent variables. As a result, for a fixed mass m the momen- 
tum distribution is Gaussian and p ~ \pm. From Eq. (10), the mass distribution is [27] 
P{m) = (m)~ 2 exp ( — m/ (m)). Direct integration of Eq. (30) shows that the momentum 
distribution is also purely exponential, P(p) = ((m)(|p|)) 1 exp(— with the typi- 
cal momentum, (\p\) = \J (m)/2. On the other hand, the velocity distribution is algebraic 
for large velocities p(v) ~ M~ 3 - 

It is interesting to compare these predictions with numerical simulations. Carnevale 
et al. [26] established a scaling behavior of Eq. (29) heuristically and confirmed it numer- 
ically. They also reported a distribution that is reminiscent of Eq. (30) . Their numeric 
form resembles Eq. (30) in that the mass distribution is exponential and in that for a 
fixed mass the velocity distribution is Gaussian. However, there is a significant difference 
between the two forms as the simulation data suggests that the velocity distribution is 
independent of mass. This observation is an intriguing one, since for a fixed mass m the 
typical velocity is mass dependent v(m) ~ m -1 / 2 . Jiang and Leyvraz [28] also studied 
the mass distribution and found that it is singular near the origin, P{m) ~ m -1 / 2 for 
m <^ (m), in contradiction with our approximate theory. Curiously, the mass-momentum 
distribution contains an identical singularity. However, this singularity disappears when 
the momentum is integrated out. 
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A similar line of reasoning applies in d dimensions. We assume that the density 
of an aggregate is constant, and thus, as the aggregation process evolves, the size of an 
aggregate grows indefinitely. Initially, only monomers with unit momentum occupy the 
system. Since the collision rate is also proportional to the surface area of an aggregate, 
the typical collision rate is wa d_1 , with the radius a ~ m l / d . The density satisfies the 
approximate Boltzmann equation 



dP{m, p) 
dt 



(» d_1 ) (J2 P « P) p (™ -m',p-p')- 2P(m, p) p')) , (31) 



where p is the <i-dimensional momentum. Repeating the above analysis yield the leading 
scaling behavior for the mass and the momentum, 

mr ^ t 2d/(d+2) and | p | ^ t d/(d+2)_ (32) 

The distribution function is a simple generalization of Eq. (31), 

P( m ,p )Mm >-^M-(^-^), < 33 > 



(m) 

with the average mass (m) ~ ^2d/(d+2)_ (32) suggests that ballistic aggregation has 
no upper critical dimension. Hence, it is not clear whether our approximation holds in 
sufficiently large dimension, as is the case for diffusive aggregation. In a recent study 
of the two-dimensional gas of sticky particles [30], some deviations from the mean-field 
predictions have been observed. The simulations [30] revealed that the growth exponent 
varies with the initial density. However, an exponential mass distribution and a Boltzmann 
energy distribution were found, in agreement with mean-field theory. We conclude that 
despite the crude nature of the approximation, it provides good estimates for the leading 
asymptotic behavior as well as the various distribution functions. 

Steady state can be achieved by adding particles to the system with rate h. We 
consider homogeneous and isotropic input of particles with unit mass and unit momentum. 
From Eq. (21), the distribution function reads, 

P( m , p ) = v^ m -(d+3)/2 exp c_^p!!y (34) 

V 2m / 

Furthermore, the mass distribution is given by N~ 1 P(m) ~ m -3 / 2 , with the density 
N ~ \fh. Note that the velocity kernel (v) is not important in the steady state since 
(v) oc Yl P(m)v(m) ~ J m~ z ' 2 m~ 1 ^ 2 is finite. One can also study the relaxation properties 
of the system. The relaxation towards the steady state becomes slower and slower as 
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the input rate vanishes. Following the analysis of Eq. (22), the corresponding relaxation 
exponent is obtained, 

r(h) ~ h~P, with (3 = (35) 

However, only in the low input rate limit, h <C 1, steady state is achieved. Indeed, the 
previous description is valid in the low coverage limit, i. e., when t <^ t c with t c ~ h -1 , 
while for t > t c the space is covered by a single "superparticle" . Since the relaxation 
exponent (3 is less than one, the time scales r and t c are well separated in the low input rate 
limit. Thus, the steady state distribution of Eq. (34) provides an intermediate asymptotics 
valid for r <C t <C t c . 

In contrast, the nature of steady state caused by a spatially localized particle input is 
a truly asymptotic one. The spatial dependence of the density distribution is of particular 
interest in this problem. We consider a spherically symmetric point-like source of ballistic 
particles with unit mass and unit momentum. Let P(m, p, r) = P(m, p, r, t = oo) be the 
steady-state radial mass-momentum density. For d > 2 the reaction is in fact irrelevant 
away from the source as will become clear below. The density satisfies the convection 
equation, d{r d ~ 1 P{m, p, r))/dr = 0, and thus, the concentration decays as r -( d ~ 1 ). So in 
the inhomogeneous ballistic problem, d = 2 is a critical dimension above which the reaction 
merely leads to the renormalization of the strength of the source. 

For d < 2, away from the source, particles have typically collided many times and 
they move with a velocity close to unity in the radial direction. Thus, collisions occur with 
rate proportional to the rms velocity v. The steady-state radial distribution P(m, p, r) 
satisfies 

_L-* ( r *-i P(m> p) ) = (va*- 1 ) (J2 P(m', p')P(m - m' , p - p') - 2P(m, p) £ P(m', p') 

(36) 

It is not difficult to verify that with the transformation r d_1 P(m,p) — > P(m, p) and 
R = J r ' r'^-^dr' — > t, the steady-state equation reduces to the time-dependent equation 
(31). Therefore when d < 2, the variable R ~ r 2 ~ d plays the role of time and from Eq. (32), 
m(r) ~ R2d/(2+d) ^ r 2d(2-d)/(2+d)_ As d — > 2, the exponent describing the mass growth 
vanishes, indicating that d — 2 is indeed the critical dimension, above which the typical 
mass far from the source is constant. At the critical dimension d = 2, logarithmic behavior 
occurs, R ~ log(r), and consequently, m(r) ~ R ~ log(r). Hence, for d < 2 clustering is 
significant and the typical mass is a growing function of the distance from the source. 

To determine the mass-momentum distribution, we tacitly impose boundary condi- 
tions similar to the initial conditions of the time dependent problem, rQ~ 1 P(m, p, ro) = 
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5(m — l)5(|p| — 1). The steady-state density far from the source reads 

d(|p| — m) 2 



P(m, p, r) ~ r - (d - 1} (m(r))" 2 m- d/2 exp 



m 



(m(r)) 



2m 



d < 2, (37). 



with (m(r)) ~ r 2d(2-d)/(2+d) for d < 2 an d ( m (r)) ~ log(r) for d = 2. For a fixed 
mass, the momentum distribution is Gaussian and as a result, the rms momentum is 
characterized by y / m(r). By integration of the mass-momentum density, one finds the 
concentration, N(r) ~ l/(r( d_1 )(m(r))). As a result, N(r) ~ r (2-5d+d 2 )/(2+d) for d < 2? 
and iV(r) ~ l/(rlog(r)) for d = 2. 

Returning to the time dependent problem, the steady state solution holds for r < t 
only, while for a larger the space is essentially empty. It is useful to estimate the total 
number of clusters N(t) ~ f* r d " 1 A^(r). For d < 2 one finds Af(t) ~ t(2d 2 -3d+2)/(d+2) ? and 
for <i = 2 one has M{t) ~ t/log(t). In the ballistic regime, d > 2, collisions do not cause 
a significant reduction in the number of clusters, and the total number of clusters grows 
linearly in time. 

To summarize, we write the leading asymptotic behaviors of the mass 



m(r) 



r 2d(2-d)/(2+d) d < 2 . 

log(r) d = 2; 

1 d> 2, 



(38) 



the density 



iV(r) 



r -(5d-d 2 -2)/(2+d) 

r _1 [log(r)] -1 

r -(d-i) 



d < 2; 
d = 2; 
d > 2, 



(39) 



and the total number of clusters 



t (2d 2 -3d+2)/(2+d) ^ < 2 . 

t/log(t) d = 2'; 

t d > 2. 



(40) 



The typical rms momentum behaves as ^Jm{r). Numerical simulations agree with the 
above in d = 1 [26]. It will be interesting to test these predictions in higher dimensions. 
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IV. Diffusive Aggregation 

In this section, we consider the diffusive driven aggregation processes involving k 
distinct species. Each of the species masses is conserved and thus there are k conservation 
laws. We apply the rate theory described in Sec. II and also investigate low dimensional 
systems where the rate equation description is expected to fail [15]. In all dimensions we 
find that in addition to the typical mass scale, there exists an additional mass scale. 

In spatial dimensions larger than the critical dimension, d c = 2, Smoluchowski rate 
theory is exact. In dimensions lower than the critical dimension, spatial correlations are 
significant asymptotically. Particles are repelling each other, and the effective reaction 
rate, k, depends on the density [20] 



N 2 / d ~ 1 d<2; 

l/|log(iV)| d = 2- (41) 
1 d> 2. 



This reaction rate can be absorbed into a suitably defined time variable T, 

With this novel time, Smoluchowski rate equation reduces to Eq. (1). 

For d < d c , this approximation yields erroneous results for the mass distribution. 
Nevertheless, it produces correct asymptotic scaling behavior for quantities such as the 
typical mass. Thus for the general case d < d c we just quote the leading asymptotic 
behaviors while in one dimension we provide exact results. Let us consider the case where 
there are two species, say A and B. A cluster is characterized by the respective masses 
of its components, tua and mg. The typical total mass is given by m A + m B ~ T, or 
equivalent ly, 

f t d / 2 d<2; 
m A + m B ~ < t/\og(t) d = 2; (43) 
U d>2. 

On the other hand, since the mass difference is a Gaussian variable, for a fixed total mass, 



one has \m A — ms| ~ \JmA + nis • This additional mass scale can be also expressed in 
terms of time, 

( t d ^ d<2; 
\m A - m B \ ~ < y/t/\og{t) d = 2; (44) 
[Vt d>2. 

Note also that the present two- mass aggregation process can be mapped onto a two-species 
aggregation-annihilation process with two conservation laws [23]. Indeed, the above results 
agree with analytical and numerical findings of Refs. [23,32]. 
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In one dimension, it is possible to derive a complete analytical solution. We consider 
a linear lattice on which point clusters hop randomly from site to nearest neighbor sites. 
We assume that the diffusion coefficient D does not depend on cluster's mass. We also 
assume that initially each site is occupied by some monomer, Aoy B with equal probability. 
The multivariate distribution function can be expressed in the form of Eq. (11), with the 
sum variable given by the Spouge's solution [19] of the ordinary one dimensional diffusion- 
controlled aggregation problem, 

P(m+,t) = e- 4Dt [/ m+ _ 1 (4 J Dt)-/ m++1 (4 J Dt)], (45) 

where I n denotes a modified Bessel function. Introducing now the scaling variables, M + = 
m + /\ / 8Dt and M_ = m_/(8Dt) 1 / 4 , one can rewrite the solution in a convenient scaling 
form 

P(mi, to 2 , t) ~ t" 5/4 $(M+, M_), (46) 

with the scaling function <E> 

2 

$(x,y) = v^exp(-a: 2 - |^). (47) 

Finally, we consider diffusive aggregation in a system with a steady spatially localized 
monomer input. We assume that monomers of all types are added at random with an equal 
rate. Making use of exact and scaling results for the corresponding ordinary aggregation 
[20], we solve for our case. This inhomogeneous system is characterized by two critical 
dimensions, the usual "homogeneous" critical dimension d c = 2 and an additional critical 
dimension d c = 4 which demarcates the pure diffusion regime d > 4 (particles do not 
affect each other far away from the source) and the diffusion-reaction regime d < 4. The 
system approaches steady state as t — > oo. In the diffusion- reaction regime, the density 
distribution function approaches a power-law in r, where r is the distance from the source. 
The borderline cases d = d c = 2 and d = d c = 4 should be treated more carefully since 
logarithmic factors appear. We write the final results for the typical total mass 



< r 2 d<2; 
r 2 / log(r) d = 2; 



rriA + tub 



r 



4- 



2 < d < 4; (48) 



log(r) d = 4; 
1 d > 4. 



For a fixed total mass tua + ms, the mass difference is again Gaussian and consequently, 
\mA - m B \ ~ V m A + m B - 
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V. Conclusions 



In summary, we have investigated irreversible aggregation with many conservation 
laws. The solution to this process is characterized by the Gaussian statistics of the fluc- 
tuations in a given conserved quantity. The process is thus governed asymptotically by 
two size scales. Application to a "sticky gas" suggests a Boltzmann velocity distribution 
for a fixed mass. In addition, the mass distribution is exponential. By comparing our 
predictions with available numerical results we have found that the present approximate 
theory gives a good description of the sticky gas. We have also investigated steady-state 
properties by introducing a localized source. We have observed two different behaviors, 
the ballistic- reaction regime for d < 2, and the pure ballistic regime for d > 2. Application 
to diffusive aggregation with more than one species also exhibit a novel "diffusive" size 
scale. 

This study suggests different avenues for further investigation. The theory might be 
applicable to problems such as catalysis and chemical reactions with many species. It 
will be also interesting to analyze the rate equations with more realistic reaction rates. An 
important question to be addressed is how robust is the Gaussian nature of the fluctuations 
statistics. It is plausible that spatial correlations introduce nontrivial internal arrangement 
of the clusters, leading to more complicated statistics. It is also plausible that even on 
the rate equation level but with the reaction rate not expressible as a function of the sum- 
variables only, different asymptotic behavior emerges. This could explain why for a dual 
fragmentation process the multivariate generalization produces an infinite set of scales [33- 
35] compared to the two scales in the models of multivariate aggregation we have examined 
in this study. 
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